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ABSTRACT 

We present a set of four Gemini-North GMOS/IFU observations of the central dis- 
turbed regions of the dwarf irregular starburst galaxy NGC 1569, surrounding the 
O ■ well-known super star clusters A and B. This continues on directly from a compan- 

ion paper (Westmoquette et al. 2007a), in which we describe the data reduction and 
^ ■ analysis techniques employed and present the analysis of one of the IFU pointings. By 

' decomposing the emission line profiles across the IFU fields, we map out the prop- 

erties of each individual component identified and identify a number of relationships 
and correlations that allow us to investigate in detail the state of the ionized ISM. Our 
^ ' observations support and expand on the main findings from the analysis of the first 

IFU position, where we conclude that a broad (^400 kms ^) component underlying 
the bright nebular emission lines is produced in a turbulent mixing layer on the surface 
of cool gas knots, set up by the impact of the fast-fiowing cluster winds. We discuss 
the kinematic, electron density and excitation maps of each region in detail and com- 
OO ■ pare our results to previous studies. Our analysis reveals a very complex environment 

\ with many overlapping and superimposed components, including dissolving gas knots, 

rapidly expanding shocked shells and embedded ionizing sources, but no evidence for 
organised bulk motions. We conclude that the four IFU positions presented here lie 
well within the starburst region where energy is injected, and, from the lack of sub- 
stantial ordered gas fiows, within the quasi-hydrostatic zone of the wind interior to the 
sonic point. The net outfiow occurs at radii beyond 100-200 pc, but our data imply 
^ . that mass-loading of the hot ISM is active even at the roots of the wind. 

Key words: galaxies: evolution - galaxies: individual: NGC 1569 - galaxies: ISM - 
galaxies: starburst - ISM: kinematics and dynamics. 
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1 INTRODUCTION 

Understanding feedback effects between massive stars, star 
clusters and the interstellar medium (ISM) is fundamental to 
the study of galaxy evolution. This is the second paper pre- 
senting integral field unit (IFU) observations of ionized gas 
in the dwarf starburst galaxy NGC 1569. Here we focus on 
the state of the ISM in the centre of the galaxy, by studying 
in detail four regions showing characteristics of wind-clump 
interaction. 

The dwarf galaxy NGC 1569, is known to have un- 
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dergone a significant recent global starburst event of 
^,100 Myr duration, the peak of which ended 5-10 Myr ago 
(|Greggio et al.|[r998l l . During this event the galaxy witnessed 
at least 10^ supernovae (SN) explo sions, and the formation 
of many star clusters (|Anders et al.li2004il , inclu ding the two 
well-known supe r- star clusters (SSCs) A and B (|Ableslll97ll : 
lAro k Sandlg3ll985l ). The result of the energy input into 
the ambient medium can clearly be seen: NGC 1569 has a 
highly turbul ent ISM, e xhibiting ionized b ubbles, arcs and 
filaments ( Hunter et al.| [l993: H eckman et a l. 1995), and ex- 
hibits structure down to the resolution limit of every obser- 
vation made so far. 

HST images, embodying the highest resolution obser- 
vations currently available, illustra te well the complexity of 
the small-scale ISM in NGC 1569. iBuckalew k Kobulnickvl 
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(|2006l ) used HS T /WFFC2 narrow-band imaging to examine 
the detailed structure of the ISM and the underlying ionizing 
mechanisms. They used the ratios of [S ii]/Ha and [O III] /H/3 
and t he theoretical "maximum starburst line" ([Kewlev et alJ 
l200lh to determine a threshold between photoionized and 
non-pho toionized (assumed to be shock excited) emission 
(see also ICalzetti et al"]|2004h . Identifying pixels in this way, 
they found structures of non-photoioinized points across the 
face of the whole galaxy. As well as fin ding a great deal of 
activity surrounding H II region No. 2 (| Waller! 199 ih . they 
found "walls" or arcs of non-photoionized points near SSCs 
A and B, indicating interaction between their respective 
cluster winds and the surrounding ISM. Another cluster of 
points was found between SSC A and cluster 10, indicating 
a further possible wind-wind interaction between these two 
clusters. 

The distribution of He ii A4686 emissio n in NGC 1569 
was r eported in an earlier contribution by iBuckalew et al.l 
(|200Q| 1. where its signature was used to search for Wolf- 
Ray et stars. Since Hell emission can also be produced in 
fast s hocks (|Garnett et all Il99lh , iBuckalew Kobulnickvl 
(|2006l l conclude that where non-photoionized points were 
found to be coincident with Hell emission, this is strongly 
suggestive of shock- ionized gas. In a starburst environment, 
shock-excited emission dominates when the level of pho- 
toionization from very young stars starts to decline, and 
the energy input from stellar winds and SNe (which drive 
the shocks) is constant. This co rresponds to the period 
6-40 Myr (iLeitherer et al.1 Il999l ) and is consistent with 
the post-starburst classifica tion and the typical ages of 
star clusters in NGC 1569 jG onzalez- Dekado et al.1 Il997l : 
lAnders et alll2004l : lAngeretti et al. 20Q5h . However, pock- 
ets of on-going star formation must still persist in smaller 
areas of the galaxy (e.g. H il region 2) to provide the contin- 
ue d photoioni zed component. This situation was confirmed 
bv lTokura eTal.1 ()20oi ). who, using observations made with 
the Subaru telescope, detected a number of strong mid-IR 
sources in NGC 1569, the most prominent of which (MIRl) 
is coincident with Hii region 2. This and other detections 
indicate continued, deeply embedded star formation in H ii 
region 2 and the surrounding regions. 

As well as expanding bubbles and shells of shocked gas, 
evidence of the past starburst a ctivity can a l so be found 
through other, more direct ways. ICreve et al ] (|2002') iden- 
tify a number of radio supernovae (SNe) and supernova rem- 
nants (SNRs) distributed throughout the NGC 1569 disc 
using MERLIN and VLA 1.4 GHz and 5 GHz observations. 
The number of non-thermal sources identified throughout 
the central regions indicate that a great deal of star forma- 
tion has taken place in the SSC A/B region in the past. 

The collective effect of the massive stars in the clusters 
formed in the starburst has injected a great deal of kinetic 
energy and momentum into the surrounding gas, together 
with chemically processed material and copious ionising ra- 
diation. This has resulted in an impressive gal axy- wide out - 
fiow seen to extend over kpc scales in both Ha (iMartinll 19981 ) 
and X-rays (Heckman et al.l ll995l : iMartin et al.ll2002l ). 

iHe ckman et al. (1995) also presented deep Ha long- slit 
spectroscopy of the disc and filament system, and found the 
Ha line-profiles within the galaxy core to be "distinctly non- 
Gaussian" . Decomposing them using mult i- Gaussian profile 
fitting allowed them to identify two kinematically distinct 



subsystems: a bright and centrally located 'quiescent sys- 
tem' centred on SSC A, with line widths of FWHM = 30- 
90 kms~^ and radial velocities within 30 kms~^ of their 
adopted Vsys (—68 kms~^); and a fainter and extended 'high- 
velocity system' with FWHMs ^ 150 kms~^ and velocities 
up to 200 kms~^ relati v e to t heir Vsys- This confirms the 
finding of iTomita et al ] (Il994l ) who also report two simi- 
lar velocity systems. Furthermore, at the location of SSC 
A, they found the Ha profile to have weak but very broad 
wings. The investigation of this enigmatic, underlying broad 
component formed one of the main focusses of the first part 
of th is study described in Westmoquette et al. (2007a, here- 
after IPaper j). 

In IPaper j . we focussed on one IFU position, covering 
the area surrounding cluster 10 just to the east of H ii region 
2. Through accurately modelling the observed line profiles 
we were able to map the emission line properties of each line 
component across the IFU field-of-view (FoV). We found 
that the observed nebular emission lines are composed of 
a bright, narrow component and an underlying, very broad 
(<300 kms~^) component. By mapping out, for the first 
time, the properties of both the narrow component (which 
we refer to as CI) and the broad component (referred to as 
C2) in two dimensions, we were able to perform a detailed 
investigation of each component and their relationship to 
one another, and develop a new semi-unified model of the 
state of the ionized gas clumps in the ISM to explain the 
origin of the line profiles. We concluded that the bright, nar- 
row component (CI) represents the ambient turbulent ISM, 
where its motions (line width) result from a convolution of 
the general stirring effects of the starburst and gravitational 
virial motions. We found that the width of C2 is both highly 
correlated with the intensity of CI and anti-correlated with 
its own fiux. This led us to conclude that the most likely 
explanation f or C2 is that it originate s from a turbulent 
mixmg laver (iBegelman Fabianlll990l : ISlavin et al.1ll993l : 
iBinette et al.lll999l ) set up on the surface of the gas clumps 
by the effects of the strong cluster winds. In this interpreta- 
tion, the turbulent gas is then stripped from this interface 
layer through hydrodynamical effects such as thermal evap- 
oration and mechanical ablation as it becomes entrained in 
the wind fiow, resulting in a diffuse, turbulent gas compo- 
nent pervading the whole region. 

In this paper, we extend our investigation to the re- 
maining three IFU fields. In Section [2] we describe the ob- 
servations taken and summari se the data-reduction process 
described in detail in IPaper j . Then, in Section [3l we per- 
form a detailed study of each remainin g IFU position and, 
armed with our findings from IPaper j , investigate how our 
prior conclusions apply to the other regions within the star- 
burst core. 

Previous studies of the ISM in starburst galaxies and 
their resulting outfiows have tended to concentrate on large- 
scale features (mainly due to instrumental restrictions). 
Hence there exists a paucity of data of sufficient quality 
and resolution to compare directly to the latest generation 
of high-resol ution hydro dynamic starburst-d riven outfiow 
mode ls (e.g. IStevens fc^H artwell 2003; Te nor io-Tagle et aP 
I2OO3I V By examining the collective state of the ISM as shown 
by all four fields together in Section 2] we attempt to begin 
filling in these observational deficiencies. Finally our results 
and conclusions are summarised in Section [S] 
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Table 1. Gemini GMOS observations 



Pos. 


Date 


Coordinates 


Exp Time 


No. 




(J2000) 


(s) 


1 


7/10/04 


04^30^48^38 64° 51^06(^5 


1600 X 4 


2 


7/11/04 


04^30^4^827 64°51'03'/2 


1680 X 4 


3 


7/10/04 


04^30^48!04 64° 50^54(^8 


940 X 4 


4 


8/11/04 


04^30^49^58 64°50'58('9 


1680 X 4 



In a third paper (Westmoquette et al. in prep.; Paper 
III) we will examine the properties of the galactic outflow as 
a whole using deep Ha imaging and IFU observations and 
attempt to create a unified physical picture of the outflow 
state by connecting what we find here in the central regions 
to the outer halo. 



2 OBSERVATIONS AND DATA REDUCTION 

In October and November 2004 we obtained Gemini-North 
Mult i- Object Spectrograph (GMOS) data covering four re- 
gions i n the nucleus of NGC 1569 wi th the integral field unit 
(IFU; lAllington-Smith et al.1 l2002l ) (programme ID: GN- 
2004B-Q-33, PI: L.J. Smith). A f ull descri ption of the in- 
strument and the set-up is given in IPaper j . but briefly, the 
IFU in one- slit mode gives a field-of-view of 5 x 3.5 arcsecs 
(which corresponds to approximately 50 x 35 pc at the dis- 
tance of NGC 1569) sampled by 500 hexagonal contiguous 
fibres of ^'!2 diameter. A separate offset block of sky fibres 
provides a simultaneous sky view. We took four exposures at 
each position with integration times between 940 and 1680s 
each, and used the R831 grating to give a spectral coverage 
of 4740-6860 A and a dispersion of 0.34 Apix-^ Table □ 
Usts the coordinates and exposure times for each position, 
and their respective locations are shown in Fig. [1] 

Position 1 was chosen to intercept the (possibly) high 
velocity outfiow emanating from the central regions of the 
galaxy. Position 2 is coincident with cluster 10 ([Hunter et al.l 
I2OOQI I and covers some of the Ha gas on the easternmost 
part of the ionized complex centred on Hii region 2 (IWallerl 
Il99lh . A detailed description and analysis of this position, 
including a discussion of the relat ionship between cluster 
10 and its environment, is given in IPaper 11 . Position 3 was 
placed to the south of SSC A, and covers a number of com- 
pact ionized knots. VLA-8, a non- thermal radio source which 
iGreve et al.l (|2002h attribute to being low surface brightness 
SNR, is located at the very top of the field-of-view. IFU 
position 4 was chosen to focus on an ionized knot to the 
east of SSC A, north of SSC B, and south of cluster 30. We 
suspected this region might show signs of the interaction of 
winds from these three star clusters . This knot a lso contains 
the non-thermal radio source M-5 (iGreve et al. l l2002 ). 

A full description of the data reduction steps is given 
in IPaper j . Here we briefiy summarise the procedure and 
describe the changes in the differential atmospheric refrac- 
tion (DAR) correction we had to make for the remaining 
fields. Basic reduction was performed following the standard 
Gemini reduction pipeline (implemented in irafQ). First, a 

^ The Image Reduction and Analysis Facility (iraf) is dis- 



Table 2. Differential atmospheric refraction (DAR) corrections 
(YijS to Ha). 1 spaxel has a diameter of 0^18. 



Position 


Spaxel Shift 


1 


1.50 east, 2.50 north 


2 


0.39 west, 2.98 northt 


3 


0.50 west, 2.50 north 


4 


1.15 west, 2.75 north 



measured using the ce ntroid fit ting routine described in 
iPapeT^ . 

trace of the position of each spectrum on the CCD was pro- 
duced from the fiat field. Throughput correction functions 
and wavelength calibration solutions were created and ap- 
plied to the science data, before the individual spectra were 
extracted to produce a data file containing 750 reduced spec- 
tra. The final steps were to clean the cosmic-rays, subtract 
an averaged sky spectrum (computed from the separate sky 
field) , and apply a fiux calibration derived from observations 
of the standard star G191-B2B. Separation of the sky spec- 
tra from the science data resulted in a data file formed of 
500 fully reduced spectra, with a resolution varying between 
FWHM = 74±5 kms~^ at the blue end (with an average 
S/N of 4.6 in the continuum), to 59 ±2 kms~^ at the red 
end (with an average S/N of 14.5). 

DAR is the effect of varying light refraction by the 
Earth's atmosphere at different wavelengths. In order to 
correct our data for DAR, we had to apply a spatial in- 
terpolation to resample the data from the hexagonal fibre 
arrangement into con tiguous squares. The methods are fully 
described in IPaper A together with how the DAR shift was 
measured in position 2 by tracing the position of cluster 
10 through wavelength space to produce a DAR curve. Un- 
fortunately none of the other fields contain an unresolved 
source suitable for tracing the effects of DAR, so we had to 
employ a visual inspection method t o determine th e DAR 
shift. The Euro3D Visualisation Tool (|Sanchezll20o3 ) allows 
a spatial intensity map of a specific wavelength range for 
each IFU field to be plotted. After selecting two maps of the 
field isolating the H/3 and Ha emission lines respectively, 
we measured the offset in position of a representative bright 
gas knot between these two known wavelengths. Assuming 
that these measurements represent the wavelength extrema 
(a reasonable approximation for our set-up), we normalised 
the curve output from the automatic DAR tracking algo- 
rithm to our measured values in order to get the DAR cor- 
rections for these positions. The calculated offsets from H/3 
to Ha for each field are given in Table [2] 



3 PROPERTIES OF THE IONIZED GAS 

3.1 Decomposing the line profiles 

In order to quantify the gas properties observed in each IFU 
field, we fitted multiple Gaussian profile models to each 

tributed by the National Optical Astronomy Observatories which 
is operated by the Association of Universities for Research in 
Astronomy, Inc. under cooperative agreement with the National 
Science Foundation. 
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Figure 1. HST F656N image of t he centre of NGC 1569 showing the positions of th e four IFU fields. A number of the most 
prominent star clusters (|Hunter et al.|[200oh and radio continuum sources (iGreve et al.l 120021) have been labelled in black and white, 
respectively 



emission line us in^ a progra m written in IDL called pan 
(Peak ANalvsis: lDimeoll2005l V A detailed descripti on of the 
program and our modifications is given in lPaper 1 together 
with the statistical techniques we employed to determine the 
optimum number of Gaussian components to fit to each line. 

The general shape of the emission line profiles observed 
in all four positions is a convolution of a bright, narrow com- 
ponent overlying a faint, broad component. In£aper I, we 
describe our methodology and nomenclature for defining ini- 
tial guess parameters for single-, double- and triple-Gaussian 
fits. To summarise, each line in each of the 500 spectra was 
fitted using a single-, double- and triple-Gaussian compo- 
nent initial guess, where line fluxes were constrained to be 
positive and widths to be greater than the instrumental con- 
tribution. When a double component flt was deemed most 
appropriate, flts were always specifled with the first Gaus- 
sian as the narrow component (which we refer to as CI), and 



the second Gaussian as a broader component (referred to as 
C2). Never did we find a double-peaked broad line. For triple 
component fits, an additional component (C3) was specified 
with a initial guess assigning it to a supplemental narrow 
line at the same wavelength as the main narrow line (CI). 
We employed the statistical F-test, giving the minimal in- 
crease of the fit ratio that would be required for one fit to 
be statistically distinguishable from another, together with 
a number of additional physical tests, to determine the best 
number of components to fit to each line. Clearly the quality 
of the fits depends significantly on the S/N of the emission 
line. We estimate the uncertainties on the fit properties in 
I Paper ll (section 3.1); to summarise, we found that for the 
Ha line flux, uncertainties vary between 0.5-10 per cent for 
CI (for high-low S/N lines), 8-15 per cent for C2 and 10-80 
per cent for C3. FWHM errors vary between 0.25-2 kms~^ 
for high-low S/N CI lines respectively, 4-15 kms~^ for C2 
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and 15-20 kms~^ for C3. For the radial velocity, the error 
in CI varies between 0.1-3 kms~^, 2-5 kms~^ for C2 and 
10-30 kms~^ for C3. We estimated that the addition of a 
third component (where required) increases the errors in the 
CI and C2 flux by 5-10 per cent, FWHM by --5 kms"^ and 
the central velocity by ^3 kms~^. 

Once the profile of each line in each spectrum was de- 
composed using this met hod, we u sed the visualisation tool, 
DAISY, also described in IPaper A to map the spatial vari- 
ation of the Gaussian fit properties and derived values for 
each CMOS position. The flux, width and velocity maps of 
all the emission lines look qualitatively similar, so we con- 
centrate prim arily on the Ha line since it has the highest 
S/N (see also IPaper ll ). In the following sections, we dis- 
cuss the observational characteristics of IFU positions 1, 3 
and 4 in turn with reference to the aforementioned maps. In 
each case, we begin by describing the results, including any 
trends and correlations seen, then attempt to interpret what 
we find and describe how it relates to the central region as 
a whole. 

Following IPaper A we adopt a systemic velocity, Vsys — 
—80 kms~^ throughout the analysis. 

3.2 CMOS Position 1 

Position 1, the most northerly of our four positions, coincides 
with two large knots of gas. The 'swept-back' morphology 
of these (and other nearby) knots suggests that they have 
been shaped by interactions with the material being evacu- 
ated from the region surrounding SSC A, and this position 
was specifically chosen to look for evidence of these interac- 
tions. Figs [JHl] show maps of the Gaussian properties (flux, 
FWHM and radial velocity) for each line component over 
the whole FoV. The two gas knots show up clearly in the CI 
flux map (Fig. [21 left panel), and provide convenient refer- 
ence points for interpretation of the other maps. Unlike the 
other positions, we require not only a bright narrow com- 
ponent (CI) and a broad underlying component (C2) to fit 
the Ha line profile, but a third Gaussian was needed to fit 
the integrated Une profile over a significant proportion of 
the field. Because of the additional uncertainties introduced 
by fitting a third component, we precede our discussion of 
the results with a brief description of how this has affected 
the emission-line property maps and how this affects our 
interpretation of them. 

The discontinuous region towards the central-north- 
west of the field seen in the C2 flux and FWHM maps is 
caused by the sudden change from fitting with triple- to 
dou ble- Gaus sian profiles. According to the F-test method 
see IPaper J ), the ratio in this region became no longer 
sufficient to justify a triple- Gaussian fit. The artificially high 
flux and low FWHM in C2 has resulted from fitting two com- 
ponents to an intrinsically three-component profile, since C2 
now includes some contribution from the ignored third com- 
ponent. The continuous nature of the radial velocity maps 
(which are less affected by poor S/N) proves that the com- 
ponent assignments were correctly made over the whole field 
and that we can safely ignore only the discontinuous region 
in the maps in the following discussion. 

The CI flux map (Fig. H left panel) highlights the 
two gas condensations, which, without any dynamical in- 
formation appear very similar. From the equivalent FWHM 



(Fig. [3l left panel) and radial velocity (Fig. |4l left panel) 
maps, it is clear that the two do not share the same charac- 
teristics. The southern knot stands out for having the broad- 
est (350-400 kms"^) C2 widths of all the fields (see Fig. [S] 
centre panel) and for having strong associated emission in 
C3. However the peak in C3 emission (Fig. (2] right panel), 
found just to the west of this knot, does not have a counter- 
part in the HST F656N image (Fig[T]). Between this bright 
area and the CI knot, a faint semi-circular structure can be 
identified (indicated on Fig. [2] with a white line), with a con- 
centric morphology centred on the southern CI knot. The 
FWHM of C3 increases radially towards the southern knot 
following this concentric pattern, while its radial velocity de- 
creases slightly, resulting in a velocity difference between CI 
and C3 at the knot's location of 45-55 kms~^. Fig. [5]shows 
two examples of Ha line profiles from this position, including 
one from this bright C3 region (the location of the spaxels 
from which the profiles are from are marked with circles on 
the right-hand panel of Fig. [2]). These examples clearly show 
the existence of three components; the high quality of the 
fits can be qualitatively assessed through inspection of the 
residual plots shown below the main plot in each case (see 
IPaper ll for an explanation of how the residuals are calcu- 
lated). In the southern knot, there is a very good correlation 
between CI flux and C2 FWHM (i.e. the brightest CI knots 
have the widest C2 FWHM), and anti-correlation between 
C2 flux and C2 FWHM (i.e. the brightest regions in C2 
have the narrowest C2 width), as we found for position 2 
(jPaper ll) . 

The northern CI knot has a redshifted velocity quite 
distinct from its surroundings (+5 kms~^, compared to 
the field average of ^ — 15 kms~^), but the velocity of C3 
at this point is even more redshifted (-+100 kms"^). At 
the location of the knot, the C1-C3 velocity difference is 
^-^90 kms~^ (an example of a line profile from this region is 
shown in Fig. (5] right panel). If this difference results from 
expansion, then CI represents the near side of the northern 
clump, and it must be located further along the line-of-sight 
than the southern knot and be expanding with a velocity, 
Ve^^ ~ 45 kms~^ (half the C1-C3 velocity difference). 

The properties of the inter-clump gas are generally 
fairly uniform. The only unusual relationship is that CI has 
a blueshifted (^ — 15 kms~^) mean velocity, where as C2 an d 
C3 are both redshifted relative to Vsys (-80 kms~^ IPaper ll ). 
A region with distinctly different characteristics is found in 
the far south-west. Here the gas becomes rapidly redshifted 
in both CI and C3, and is coincident with an area of very 
broad ('^100 kms~^) CI emission. This corner is faint in all 
three line components, and by comparing its location to the 
HST image, appears to be part of the low surface-brightness 
cavity surrounding SSC A. 

The electron density of the ionized gas can be derived 
from the flux ratio of [S ii]A6713/A6731, but the S/N of the 
[S ii] lines in this field were found to be too low to make any 
confident profile fits. We have therefore spatially binned the 
datacube by 2x2 in order to increase the S/N per spaxel, 
and fitted the resultant [S ii] profiles with pan as we did 
for Ha. The resulting map for CI only is shown in Fig. [6] 
(left panel). The CI electron densities across the field are 
in general below or at the low density limit (indicated by 
the single contour line) , with the notable exception of three 
regions where the densities reach a few hundred cm~^. One 
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Figure 2. Pos 1. Left: Flux map in Ho; CI (range 112-2213); centre: flux map in Ha C2 (range 4-275); right: flux map in Yia C3 (range 
4.3-391). Non-detections are represented as hatched spaxels, the x and y scales are in arcseconds offset from the lower-left spaxel, and a 
scale bar is given for each plot with units 10"-*^^ erg s"-*^ cm~^ spaxel"^. North is up and east is left. The spaxels from which the example 
Ha proflles shown in Fig. [5]were extracted are marked with circles, and the semi-circular feature discussed in the text is outlined with a 
white line. 
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Figure 3. Pos 1. Left: FWHM map in Ha CI (range 13-95); centre: FWHM map in Ha C2 (range 120-370); right: FWHM map in Ha 
C3 (range 4-135). A scale bar is given for each plot in kms"-*^, and the measurements are corrected for instrumental broadening. North 
is up and east is left. 



is found in the south of the field extending up towards the 
north-west, coincident with the bright CI and C3 emission 
associated with the southern knot; the other two are located 
in the centre and north-west of the field, apparently not 
associated with any bright emission at all. This result shows 
that the compact nature of the two knots seen in the flux 



map of CI is in fact misleading, and that the density of the 
gas in these condensations (particularly the northern knot) 
is actually quite low. 

Line flux ratios can also be used to trace the ioniza- 
tion level of the emitting gas. In Fig. [6] we show the map of 
[S ii]A6717+A6731/Ha. Since not all individual components 
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Figure 4. Pos 1. Left: radial velocity map in Ha CI (range —24-10); centre: radial velocity map in Yia C2 (range —28-62); right: radial 
velocity map in Ha C3 (range —3-100). A scale bar is given for each plot in kms"-*^ (heliocentric, relative to the systemic velocity of the 
galaxy, Vsys = — 80 kms"-*^), where zero is marked with a line. North is up and east is left. 



Spaxel No = 33. 




Spaxel No = 331. 





Figure 5. Pos 1. Examples of Ho; fits for spaxels 33 and 331 (positions marked with circles on Fig. [2])- These examples illustrate the 
presence of three significant line components and their contribution to the integrated line profile. 



could be identified in all lines, we calculated the ratios from 
the sum of the flux in all observed line components. Fur- 
thermore, this diagnostic is less sensitive to low S/N since 
we sum rather than divide the flux from the two [S ii] lines, 
and we have plotted this map without having to bin the dat- 
acube. The lowest ratios are found towards the north-east 
of the field, corresponding to the northern Ha knot, and in 
the south-west, roughly coincident with the peak in Ha C3 
emission. This ratio, when combined with [O iii]A5007/H/3 
(not shown here because of the low S/N of the H/3 Une) 
can be used to look for indic ations of shocks. Following the 
method described in lPaper^ for position 2, we define points 
in [O lll]/H/3-[S ll]/Ha space as containing a contribution 
from non-photoionization processes when they fall above the 
'maximum starburst line' (|Kewlev et al.ll200lh. All points in 
this field lie below this line. We noted in IPaper ll that the 



low metallicity of NGC 1569 may mean that this threshold 
for non-photoionized (shocked) emission is in reality lower 
than the generalised 'maximum starburst line'. We therefore 
conclude that emission in this field is close to being fully 
consistent with pure photoionization. A contribution due to 
shocks may be present, but is currently unquantifiable. 

3. 2. 1 Interpretation 

This position covers two dynamically different gas knots that 
form part of the clumped medium around which the central 
cavity is being evacuated. The southern half of this field 
is similar to position 2 in that we see very broad (up to 
350Jans~^) C2 widths coincident with bright CI emission. 
In I Pap er^ we discussed in detail what possible mechanisms 
might cause such broad emission coincident with brightly 
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Figure 6. Pos 1 Left: Electron density map derived from the ratio of [S ll] A6717 /A6731 for CI [a scale bar is given in units of cm~^; 
a single contour line indicates the low density limit (100 cm~^: [QsterbrodjllQSOh : right: ratio of [S ll] A6717+A6731/Hck; in logarithmic 
scale. 



emitting knots, and concluded that the most likely expla- 
nation is a turbulent mixing layer on the surface of the gas 
cloud set up by the shear forces created by the fast-flowing 
cluster winds. Thermal evaporation and/or mechanical ab- 
lation of clump material result in mass being stripped off 
the surface of the clouds and thus a highly turbulent veloc- 
ity field pervades the whole area. We note that in every case 
(here and in position 2), the broad lines are single peaked, 
meaning we are not likely to be seeing line broadening due 
to projection through a bipolar outflow. Since the data for 
this southern knot are fully consistent with the turbulence 
explanation, we can extend our conclusions taking into ac- 
count a number of other observations. 

Firstly, in general the gas is found to have a low den- 
sity, possibly explaining why the average intensity of CI and 
C2 are respectively five and three times lower than in po- 
sition 2. These low intensities have resulted in more uncer- 
tain fits and much noisier maps. A density enhancement of 
^500 cm~^ is observed in the centre of the southern clump, 
but the gas immediately surrounding it is of very low density 
(<100 cm~^). Secondly, the width of the broad component 
over the whole field is on average 50 kms~^ larger than in 
position 2. We interpret these two facts to mean that this 
southern knot is undergoing a strong interaction with the 
cluster winds from the central region (particularly SSC A) 
and is at a more advanced stage in its dissolution than the 
clumps observed in position 2. The low density of the outer 
parts of the gas knot may mean that the wind can penetrate 
more of the clump and material is being stripped off more 
easily. 

We do not see evidence for the swept-back tail morphol- 
ogy seen in the HST Ha image, but this may again be due to 
the very low density, diffuse nature of the gas or the limited 
spatial-resolution of our IFU observations. Furthermore, as 
shown by the random distribution of radial velocities mea- 
sured in this region, we cannot be sure which cluster or 



clusters are dominating the energetics in this region, so it is 
quite possible that what might appear to be a swept-back, 
bow- shaped interaction in the HST image is just part of the 
ambient, clumped gas distribution. 

Observing such a strong additional narrow component 
(C3) over the majority of the field is an indicator of the 
disturbed nature of this part of the inner 200 pc, however 
we defer a discussion of this component until Sect ion |4] where 
we consider it in context with the whole central region. 

3.3 GMOS Position 3 

Position 3 is located just to the south of SSC A, in what ap- 
pears to be a region of turbulent filamentary gas. The radio 
source VLA-8 is located near the north o f the field at coor- 
dinates a = 4^ 30™ 48!0, 6 = 64° 50' 56'.'3 (iGreve et al.ll200l 
see Fig. [1]), and is thought to be an extended {2^' = 20 pc), 
low surface brightness SNR. To the south and west, just be- 
yond the IFU FoV, there is a large, prominent, semi-circular 
shell- li ke structure outlined by a string of non-photoioinized 
points (|Buckalew Kobulnickvl [20061 ) . indicating that this 
could be part of a fast, expanding bubble. Due to its close 
projected distance from SSC A, this region looks like it may 
bear the brunt of the impacting winds from the central star 
clusters. 

The CI intensity map (Fig.[7l left panel) shows a bright 
knot in the north-west corner, with a diffuse, extended 
area of emission spreading out towards the south-east. This 
diffuse emission corresponds to part of a highly clumped, 
ragged structure seen on the HST image of Fig. [1] which 
is unresolved in our IFU data. In the C2 flux map, the CI 
velocity map and in both the CI and C2 FWHM maps, this 
area is no longer smooth and diffuse, but rather surprisingly 
contains a compact region (O'.'S x l'.'4 = 9 x 16 pc) with 
very different properties to the rest of the field. The FWHM 
of both CI and C2 (Fig. [8]) peak here at 70 kms"^ and 
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Figure 10. Pos 3. Left: contours of CI FWHM (12 equally 
spaced levels from 4 to 73 kms"-*^) overlaid on the CI flux map; 
right: contours of CI radial velocity (12 equally spaced levels from 
—7 to +21 kms"-*^) overlaid on the CI flux map. Solid contours 
represent levels above the mean of the range in each case, and 
dashed contours represent levels below the mean. These two maps 
illustrate the offset of the peak of broad and redshifted emission 
with respect to the peak in the flux distribution. 



300 kms~^, respectively (corrected for instrumental broad- 
ening), both corresponding to a strong redshifted peak in 
CI (^^25 kms~"^ with respect to v^ys = —80 kms~"^; Fig. [71 
left panel). In Fig.[TOl we overplot the 01 FWHM and radial 
velocity maps (contours) on the CI flux map (greyscale) to 
illustrate how this compact region, so clearly defined in the 
line dynamics maps, only roughly corresponds to the north- 
west of the diffuse emission mentioned above, and certainly 
not with the bright knot in the far corner. The rapid peaks 
in line- width and the contrasting nature of this knot indicate 
that this could be the position of the SNR VLA-8. It is of 
approximately the right extent, but since its centre is offset 
by ^2^^ (based o n HS T coordinates) from the location given 
bv lGreve et aD (l2002l ). we cannot confirm this scenario. 

A faint diagonal strip in the C2 flux map runs from the 
north-east to south-west of the field in the opposite direc- 
tion to the CI morphology, wrapping around the unusual 
knot to the north and south, and is highly correlated with 
the distribution of narrowest C2 line-widths. This pattern 
does not have a counterpart on the HST image of Fig. [1] 
The north-eastern corner of this field is coincident with part 
of the evacuated region around SSC A (again see Fig. [1]), 
and contains the narrowest CI widths (as low as 5 kms~^; 
consistent with just thermal broadening) seen in any of the 
four IFU fields. 

The gas in the south of the field, fairly unremarkable 
in the flux or FWHM maps, stands out in the radial ve- 
locity maps of both CI and C2 as being very blueshifted 
('-^ 10-20 kms~"^ with respect to Vsys] representing the most 
blueshifted C2 in all fields). A representative profile (to- 
gether with the best-fitting model) from this region is shown 
in Fig. [11] When compared to the broader context as seen in 
the HST image, the blueshifted gas corresponds to an area 
just interior to the shell-like feature just off the south-west of 
the IFU FoV (described above), possibly indicating CI and 
C2 are both associated with the near-side of this bubble. 

Unlike positions 1 and 2, we do not see a correlation 
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Figure 11. Pos 3. Example of our flt to the Ha line in spaxel 230 
(position marked with a circle on Fig. [T]). This spaxel exemplifles 
the largest velocity difference between CI and C2 of any of the 
flelds. 

between C2 FWHM and CI flux. Instead there is a strong 
anti-correlation between C2 FWHM and C2 flux (the maps 
appear almost morphologically identical), and to a lesser 
extent between C2 FWHM and CI FWHM, although the 
respective ranges in FWHM are dramatically different. A 
generally faint CI results in this field having the highest 
average flux ratio of C2/C1 of ^^30 per cent. 

Electron densities were derived from the flux ratio of 
[Sii]A6717/A6731 but fall near or below the low density 
limit (100 cm~^) across most of the field. In order to in- 
crease S/N, we have again binned the data in a 2x2 man- 
ner, and resulting electron density map for CI is shown in 
Fig. [13] (left panel). We find the electron density increases to 
a few hundred cm~^ in the north-west of the field, coincident 
with the strongest Ha emission seen in Fig. [7] The [S li]/Ha 
(map shown in Fig. [13] right panel) and [Olll]/H/3 (map 
not shown) flux ratios are consistent with pure photoioniza- 
tion in every spaxel according to the 'maximum starburst 
line' method. However, the low metallicity of NGC 1569 may 
mean that this threshold is an overestimate and that some 
contribution from shocked emission is present. Nonetheless, 
the low line ratios are a little unexpected considering the 
proximity of the SNR VLA-8. 



3. 3. 1 Interpretation 

The striking feature in position 3 is the unusual compact 
knot that dominates the CI width and velocity maps, and 
the C2 flux and width maps. Fig. 1121 shows a magnified ver- 
sion of the region surrounding this IFU position on which 
the location of the compact knot is shown by a filled ellipse. 
As we have suggested, the sudden peak in line- widths could 
indicate t hat this is the loca tion of the SNR VLA-8 cata- 
logued bv lGreve etHI (l2002l ). but since we see no evidence 
for shocked line-ratios (as would be expected if this was a 
SNR) , and the fact that it is offset from the position quoted 
for VLA-8 make this explanation less tenable. Ergo, we can 
only conclude that this feature represents an isolated, unusu- 
ally turbulent knot of faintly emitting, redshifted gas that 
happens to lie on this sight-line. 

Except for at the location of this knot, the width of 
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Figure 7. Pos 3. Le/t; Flux map in Ha CI (log range 2.5-3.5); right: flux map in Ho; C2 (log range 2.1-3.2). The x and y scales are in 
arcseconds offset from the lower-left spaxel, and a scale bar is given for each plot with units log lO"-*^^ erg s"-*^ cm~^ spaxel"-*^. North is 
up and east is left. The spaxel from which the example Ha proflle shown in Fig. II II was extracted is marked with a circle. 
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Figure 8. Pos 3. Left: FWHM map in Ho; CI (range 4-73); right: FWHM map in Ho; C2 (range 97-295). A scale bar is given for each 
plot in kms"-*^, with values corrected for instrumental broadening. North is up and east is left. 



02 does not exceed ^-^200 kms~^ over the whole field, and 
falls as low as 100 kms~^ in the north-east quadrant. Away 
from the knot, the broadest widths are not associated with 
brightly emitting gas (as found in positions 1 and 2) im- 
plying that we are not seeing a strong interaction with the 
cluster winds in this position. The average ratio of C2/C1 
for this position is the highest of all the fields, giving an- 
other indication that the physical conditions in this field are 
different to that of the other positions. 

It would seem contradictory why at this close projected 
distance from SSCs A and B (^-^50 pc) the gas has not com- 



pletely been blown away, since large swathes of material at 
similar projected distances have been completely evacuated. 
Can we explain this using the observations? We have already 
mentioned that we see little evidence for a clump-wind inter- 
action in the form of broad line- widths or shocked line-ratios. 
Both CI and C2 become gradually blueshifted with respect 
to Vsys towards the south of this field, peaking in an area co- 
incident with the interior of the shell- like feature located just 
off the bottom of th e field, if the string of non-phot oionized 
points identified by iBuckalew Kobulnick"^ (l2006l ) are m- 
cluded as part of the shell edge. Their non-photoionized 
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Figure 9. Pos 3. Left: radial velocity map in Ha CI (range —7 to +21); right: radial velocity map in Ha C2 (range —22-34). A scale 
bar is given for each plot in kms"-*^ (heliocentric) corrected for the systemic velocity of the galaxy (—80 kms"-*^), where zero is marked 
with a line. 




Figure 12. HST F656N image detail showing the region 
around IFU position 3. The black hatch ed areas encompass non- 
photo ionized points found in the study of B uckalew Kobulnickvl 
(|2006l V and the filled ellipse marks the location of the unusual 
compact knot. 



points are shown in Fig. [12] as hatched areas, and indicate 
that the eastern edge of this shell, although not Ha bright, is 
shock-excited. This whole region, with the exception of the 
unusual isolated knot, may therefore be located at a different 
line-of-sight to the surrounding gas and be shielded to a cer- 
tain extent from the effects of the strong cluster winds. This 
would explain the relatively quiescent line- widths, and was 



a result pred icted for the turbulent mixing/ablation model 
described in IPaper A for regions not subject to the direct 
impact of cluster winds. 

3.4 CMOS Position 4 

Position 4 covers a bright knot of ionized gas to the east 
of SSC A that has a morphology somewhat resembling the 
southern knot in position 1. This condensation is coincident 
with the radio continuum source M- 5 (a "non-thermal source 
with some thermal contribution" ; iGreve et al.l l2002h , and 
appears to have an oval shell-like extension to the south in 
the HST F656N image (Fig [I]). 

Before discussing the emission-line maps for this posi- 
tion fFigs [T4| - [16p . we first comment that since a third com- 
ponent is required in a small region on the east of the field, 
its presence has resulted in high uncertainties on the fits t o 
C2 in this same region (for reasons discussed in IPaperT ) . 
This has caused a spurious and erratic region in the C2 
FWHM and radial velocity maps (outlined with a solid line 
on the following maps for clarity), so we ignore this sec- 
tion when discussing the properties of C2. Fig. 1171 shows an 
example Ha line profile from this region showing the best- 
fitting Gaussian model profiles and the fit residuals. In this 
particular spaxel, a third component is clearly required to 
accurately model the integrated line shape, but it is clear 
how the associated errors on the fit to C2 would increase as 
C3 becomes fainter thus leading to this area of inconsistent 
results. 

The aforementioned bright knot shows up clearly in the 
CI flux map fFig. ll4[ left panel) as expected, but is only ev- 
ident in the corresponding FWHM and radial velocity maps 
and the C2 flux map as part of a distinct arc of gas ex- 
tending out of the FoV the north-east, and partly towards 
the south-east. This arc has a consistently low CI FWHM 
(^40 kms~^; Fig. [151 left panel), and is blueshifted by —10 
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Figure 14. Pos 4. Le/t; Flux map in Ha CI (range 204-3768); centre: flux map in Ho; C2 (range 6.4-326); right: flux map in Ha C3 
(range 4.0-20). Non-detections are represented as hatched spaxels, the x and y scales are in arcseconds offset from the lower-left spaxel, 
and a scale bar is given for each plot with units lO"-*^^ erg s"-*^ cm~^ spaxel"-*^. North is up and east is left. The spaxel from which the 
example Ho; proflle shown in Fig. llTl was extracted is marked with a circle. 



to —15 kms~^ with respect to fgys fFig. I16[ left panel), im- 
plying that the bright knot is dynamically associated with 
a larger gas filament extending out either side of the FoV. 

The gas not associated with the filament is much fainter, 
and exhibits very different dynamical characteristics. A large 
circular region in the east of the field is distinct in CI for 
having coincident peaks of very broad (up to ^-^90 kms~^; 
representing the broadest CI of all fields) and redshifted 
velocities (up to +20 kms~^ relative to t'sys)- Around this 
location a third component is also required to fit the line 



shape (as mentioned above). From the velocity difference 
between CI and C3, we infer expansion velocities of order 
'^^exp ~ 80 kms""*^, with the fainter component of the two 
(C3) travelling at ^150 kms~^ towards us. An absence of 
increased line- widths or a third component at the location of 
the small shell seen to the south of the knot in the HST Ha 
image indicates that this structure may not be expanding, is 
expanding only very slowly, or not a shell at all. The broad- 
est C2 lines (^250 kms~^; Fig. I15[ centre panel) are found 
in the south-east and north-west, and in both locations are 
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Figure 15. Pos 4. Le/t; FWHM map in Ha CI (range 41-90); centre: FWHM map in Ho; C2 (range 60-270); right: FWHM map 
in Ha C3 (range 4-130). The region outlined with a sohd hne in the C2 map indicates where the measurement errors are high due to 
the addition of a third Gaussian to the fit. A scale bar is given for each plot in units of kms"-*^, with values corrected for instrumental 
broadening. North is up and east is left. 
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Figure 16. Pos 4. Left: radial velocity map in Hck; CI (range —14-20); centre: radial velocity map in Ho; C2 (range —9-124); right: 
radial velocity map in Ha C3 (range —156 to —86). Again, the region outlined with a solid line in the C2 map indicates where the 
measurement errors are high due to the addition of a third Gaussian to the fit. A scale bar is given for each plot in kms"-*^ (heliocentric) 
corrected for the systemic velocity of the galaxy (—80 kms"-*^), where zero is marked with a line in the left-hand panel and +50 kms"-*^ 
is marked with a dashed line in the central panel. 



14 M.S. Westmoquette et al. 



Spaxcl No = 



1^-° - f. 



50 - 









6554 


6556 


6558 6560 6562 
Wavelength (A) 


6564 65 


66 6568 


3 












1 

-1 

-2 













6554 6556 6558 6560 6562 6564 6566 6568 

Wavelength (A) 



Figure 17. Pos 4. Example of our fit to the Ha line in spaxel 338 
(position marked with a circle on Fig. ll4|) . This example was cho- 
sen to illustrate the presence of three significant line components 
found in the area on the east of the field. The third component, 
C3, is the small profile in the blue wing. 



associated with redshifted gas in CI (+5-20 kms-^). We 
mention this because these are sections of the IFU FoV that 
appear to be coincident with part of the evacuated cavity 
surrounding SSC A, and will be discussed in Section [4. II 

This field has a number of unexpected relationships 
between the components compared to the other positions. 
Firstly, this is the only position where the flux distribution 
of C2 (Fig. [Ml centre panel) bears any resemblance to the 
corresponding CI map and to the HST image. We show be- 
low how this can be explained by the unusually high density 
of the gas emitting C2. Secondly, the correlation between 
C2 FWHM and CI flux is in the opposite sense to what we 
have seen in positions 1 and 2. Here the broadest C2 lines are 
coincident with the faintest CI emission in the north-west, 
south-west and centre-east of the field. 

This is the only position for which we have been able 
to make confident double-Gaussian fits to some of the [S ii] 
lines, and derive the electron density for both components 
(mapped in Fig. [TS]) . The density of the gas emitting CI 
(left panel) is low, with only parts of the western edge co- 
incident with the bright knot rising above the low density 
limit (indicated by the single contour). We are able to de- 
rive the density for C2 (centre panel) in a diagonal band 
stretching from the centre- west to the north-east, coincident 
with the band of bright C2 emission and blueshifted CI gas 
in the north of the field. The density of C2 is on average 
much higher than that of CI, reaching up to ^600 cm~^ 
in the north-east. Since the emission intensity of recombi- 
nation Unes is proportional to the gas density-squared, the 
unusually high densities in C2 may may explain why in this 
field the C2 emission morphology bears a close resemblance 
to the corresponding CI morphology. 

Also shown in Fig. [18] is the flux ratio of 
[Sii](A6717+A6731)/Ha (right panel). The lowest ratios 
are found in the cent re- west of the field, coincident with 
the bright Ha knot and the radio source M-5. Using the 
[ O III] /H/3- [S ii] /Ha nebular diagnostic technique described 
in lPaperJ and above, we find a significant number of spaxels 
over the whole eastern half of the field contain a possible con- 
tribution from non-photoionization processes (presumably 




Figure 19. HST F656N image detail showing the region around 
IFU position 4. The h atched areas are non-photoioni zed points 
found in the study of iBuckalew Kobulnickvl (|2006l ), and the 
arrow indicates the direction of SSC A. 

shock excitation; see lCalzetti et al.ll2004f ). However, the de- 
rived low flux ratios at the location of the bright knot do not 
support a SNR origin for M-5 unless it is highly embedded. 

3.4-1 Interpretation 

Evident in both the flux maps of CI and C2, a band of 
material curves round from the north-east of the field to 
the south, including the bright knot near the western edge. 
In CI, this band has a consistently low FWHM and is 
blueshifted by 10-15 kms~^ with respect to v^ys- In com- 
parison, the ionized gas to the east of the field is signifi- 
cantly redshifted in CI, has a high- velocity ('-^150 kms~^) 
blueshifted component, and a widespread non-photoionized 
contribution to the excitation. 

Careful inspection of the HST image of Fig. [1] shows 
that this band may be the westernmost part of a frag- 
mented shell extending out of the FoV to the east. Fig. [T9l 
shows a blow-up of the HST F656N image around position 4. 
The hatched areas represent the regio ns of non-ph otoionized 
points found by Buckalew & Kobul nickvl (|2006l V illustrat- 
ing how the edges of this shell appear to be traced out 
by non-photoionized (shocked) emission. How can we re- 
late this to our more wide-spread distribution of non- 
photoion ized gas shown in Fig. [TSf The methods used 
by Buck alew Kobulnickvl for the determination of their 
points are inherently luminosity weighted, in the sense that 
only where they actually detected emission at sufficient S/N 
could they measure the line ratios. The fact that our data 
are significantly deeper, combined with the advantages of 
actually being able to fit the individual emission lines (thus 
avoiding the uncertainties of continuum subtraction), has 
meant that we have been able to accurately measure the 
line ratios for much fainter gas. The distribution of non- 
photoionized gas we derive includes their points arising from 
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Figure 18. Pos 4. Le/t* Electron density map derived from the ratio of [S ll]A6717/A6731 for CI (range 0-200); centre: electron density 
map f or C2 (range 0-7 00). A scale bar is given for each plot in units of cm~^; a single contour level marks the 100 cm~^ low density 
hmit (|Osterbrocklll989l ) on both maps. Right: flux ratio of log([S ll] (A6717+A6731)/H a) CI (log range -1.3 to -0.75). Spaxels marked 
with a cross indicate non-photoionized emission according to the iKewlev et al.l (l200lh maximum starburst line. 



the brighter regions and additional points identified from the 
lower surface-brightess regions, implying that shock-excited 
emission is more widespread than previously observed. 

We can now associate the emission to the east of the 
field with originating from within the expanding shell struc- 
ture, where the near-side wall is represented by the high- 
velocity 03, and the far-side by the brighter redshifted re- 
gion seen in CI. The velocities and brightness contrast be- 
tween the two components can be explained if the far-side 
of the bubble is expanding into a dense medium, causing it 
to slow down rapidly and emit brightly. The fast velocities 
and faint Ha emission of C3, imply that the near-side is ex- 
panding relatively unimpeded, meaning that an estimate of 
the true expansion velocity of the bubble can be made by 
comparing the radial velocity of the gas in the shell edge 
(the blueshifted band) to the high- velocity C3. We estimate 
Ve^p to be ~ 140 kms~^ by taking half the velocity difference 
between the two line components. 

The width of C2 remains between 100 and 150 kms~^ 
in the blueshifted shell edge including the position of the 
bright knot. This lack of very broad lines may indicate that 
a direct wind-clump interaction is not taking place here. Re- 
call that this clump contains the radio source M-5 emitting 
a mixture of thermal and non-thermal radiation. The pres- 
ence of a thermal signature points towards an embedded 
photoionizing source (e.g. a young star cluster) being the 
origin of this emission. Emission in the thermal IR would 
support this sc enario, but un f ortun ately the high-resolution 
observations of iTokura et al.1 (|2006l ) did not cover this knot. 
There is very marginal evidence for continuum emission at 
the location of the bright knot in our IFU data, but since 
we cannot identify any discrete stellar features (e.g. a WR 
bump) we cannot prove the existence of an embedded stellar 
source. 



4 THE STATE OF THE ISM 

So far we have examined the details of the ionized gas distri- 
bution in each field individually. However, there are insights 
to be gained from looking at the collective ionized gas prop- 
erties of all four fields at once, so as to assess the general 
trends and systematic variations across the ^-^200 pc wide 
region centred on SSC A as sampled by our IFU pointings. 
In Figs[2Tl[20]and[22l we plot the Gaussian flux, FWHM and 
radial velocity of the identified line components for all four 
positions in the three possible permutations. In the follow- 
ing paragraphs, we will discuss these plots together with the 
selected average properties listed in Table (3] in the context 
of understanding the general state of the ISM. 

Plots similar to Fig. [20] have been used success- 
fully in the past as diagnosti c tools to examine line 
broa dening in giant Hii regio ns (|Munoz-Tun6n et al ] ll996l : 
Mart inez-Delgado et al. I l2007l ). These studies find inclined 
bands of points within the FWHM-flux plots to be associ- 
ated with the presence of large-scale expanding shells within 
the H II region, and use their presence to constrain the evo- 
lutionary stage of the region. We note that similar inclined 
bands are observed in our plot, but caution that unlike the 
aforementioned studies, we decompose the integrated line 
profiles rather than fitting only a single Gaussian. The vari- 
ation in width of the single fitted Gaussian is one of the 
main ideas on which their conclusions are based, making a 
direct comparison with their conclusions difficult. 

Fig. [20] shows that the trend found in position 2, 
whereby the broadest emission is the faintest, also applies 
to the entire data-set. This implies that the most turbulent 
gas does not emit strongly in Ha. Plotting all the data to- 
gether also highlights the fact that the highest velocity gas 
is also the faintest (Fig. [22]). On average, position 2 contains 
the brightest Ha emission (as would be expected from the 
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Figure 20. FWHM (corrected for instrumental broadening) vs. log(flux Ha) (arbitrary scale) for all four IFU positions. Black symbols 
represent position 1, red position 2, green position 3, and blue position 4. 
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Figure 21. FWHM (corrected for instrumental broadening) vs. heliocentric radial velocity (with respect to Vsys) for all four IFU 
positions. Symbol colours are the same as Fig. [20] 



HST Ha image, Fig. [T]), whereas position 1 is the faintest. 
However, position 1 exhibits the broadest C2 widths, with 
an even spread extending up to ^400 kms~^. As we have 
previously mentioned, this is likely to indicate that this po- 
sition includes gas that is experiencing a strong interaction 
with the cluster winds. Position 3 has the most quiescent C2 
widths, most likely resulting from the fact that the gas in 
this field is shielded to some extent from the direct influence 
of the outflowing cluster winds. 



For position 2 (|PaDer the lower limit found for the 
width of CI (referred to as FWHMo), and its similarity to 
the integrated H i velocity di spersion for the wh ole of NGC 
1569 (FWHM ^ 35 kms~^ iMihle et aLlliooBI ). led us to 
conclude that it represents a convolution between the dis- 
turbed, turbulent ISM that has been stirred up by the ef- 
fects of the intense star- format ion and supernova shocks and 
gravitational virial motions. We could not rule out an ad- 
ditional minor contribution to its width from inevitable un- 
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Figure 22. Heliocentric radial velocity (with respect to Vsys) vs. log(flux Ha) (arbitrary scale) for all four IFU positions. Symbol colours 
are the same as Fig. [20] 



Table 3. Selected average (mean) emission-line properties for the 
four IFU fields. The velocity measurements are quoted with a er- 
ror of Icr on the mean, and do not take into account uncertainties 
on the individual data points. 



IFU 


Line 


Gaussian property 


Average value 


position 


component 






all 


C2/C1 


Flux ratio 


0.12 


1 


C2/C1 


Flux ratio 


0.11 


2 


C2/C1 


Flux ratio 


0.07 


3 


C2/C1 


Flux ratio 


0.30 


4 


C2/C1 


Flux ratio 


0.13 


all 


CI 


FWHM« 


46.9 ±8.7 


all 


C2 


FWHM 


192.0 ±47.3 


all 


C3 


FWHM 


53.1 ±30.0 


all 


CI 


Radial Velocity ^ 


-1.8±6.7 


all 


C2 


Radial Velocity 


±13.4 ±17.4 


all 


C3 


Radial Velocity 


-14.6 ±24.3^ 


1 


CI 


Radial Velocity 


-16.7±6.0 


2 


CI 


Radial Velocity 


-2.5±5.5 


3 


CI 


Radial Velocity 


-3.3±5.9 


4 


CI 


Radial Velocity 


-0.1±9.3 


1 


C2 


Radial Velocity 


±19.5±21.9 


2 


C2 


Radial Velocity 


±4.4±8.7 


3 


C2 


Radial Velocity 


±4.4 ±14.4 


4 


C2 


Radial Velocity 


±25.3 ±24.7 



^ FWHMs are in units of kms"-*^, corrected for instrumental, 
but not thermal broadening. 

^ Radial velocities are in units of kms"-*^, relative to ^'sys (= 
-80 kms-i). 

^ Average without position 4 = ±24.5 ±20.6 kms"-*^ 



resolved kinematical components along the line-of-sight. By 
plotting the FWHM results for all fields together (Figs [21] 
and l2Qp , we find this consistency in CI line- width is mir- 



rored across all positions, and that they fall in a very nar- 
row band (47 ± 9 kms~^; Tabled extending over the full 
range of intensities observed. The width of CI is therefore 
not be related to its emission strength. The clear lower-limit 
to the FWHM of CI found in position 2 is also evident in the 
other fields. In positions 1 and 3, FWHMo ^ 30-35 kms~^ 
(although is less well defined for position 3), and in posi- 
tion 4, FWHMo is slightly higher at ^45 kms"\ This level 
of consist ency ther efore supports our conclusions about CI 
reached in lPaper^ . Although in places, CI is associated with 
fast motion and/or expansion, its line- width is everywhere 
equal to or greater lower limit, FWHMo, of the gas across 
the whole central regions. 

We can identify a third component in positions 1, 2 and 
4, and from Figs [H] and [20l see that its FWHM is equiva- 
lent to that of CI (within the uncertainties) in every case. 
This confirms the conclusion that where identified, C3 rep- 
resents an additional narrow component associated with a 
secondary region of ionized gas, possibly one half of a wind- 
blown bubble. In position 4 its relation to expanding shells 
is clear, but in position 1 its significance is not so obvi- 
ous. This component occupies a prominent diagonal band 
in Fig. I21[ distinct from all the other components, between 
velocity = ±100 kms"\ FWHM = 20-30 kms"\ and ve- 
locity = 30 kms"\ FWHM = 100 kms"\ implying that 
bulk motion and turbulence influence the dynamics of this 
component to differing degrees from one half of the field to 
the other. This behaviour is mirrored in Figs[20]and Fig. 1221 
where C3 of position 1 again occupies its own unique space. 
Its FWHM seems correlated with that of C2 at the posi- 
tion of the southern knot, implying that it is affected by 
the wind-clump interaction (unlike CI), however its veloc- 
ity is spatially correlated with CI only at the position of the 
northern knot. Consequently, we conclude that in position 1, 
C3 represents a mix of overlapping kinematical systems orig- 
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inating from separate mechanisms. Where it exhibits a high 
comparative redshift (in the north), it may represent the 
far-side of an expanding structure - its properties are consis- 
tent with being p art of the high- veloc ity ki nematic subsys- 
tem id entified by iTomita et alJ (I1994I ) and iHeckman et al.l 
Where its velocity is similar to CI (in the south), it 
is likely to be associated with the wind-clump in teraction , 
since its FWHM is correlated with that of C2 fsee lPaper j ). 
In summary, our C3 data reveal hidden shells that are not 
found from morphology alone. Most of these have modest 
expansion velocities (30-40 kms~^), which are enough to 
stir up the ISM but not drive a wind or produce a hot ISM. 

As listed in Table El the average radial velocity of 
C2 is offset from the corresponding CI average by only 
^^+15 kms~^ over the four fields; the largest individual dif- 
ference being in position 1 of ^35 kms~^ for reasons dis- 
cussed above. Furthermore, ignoring the high C3 velocities 
in positions 1 and 4, the total velocity spread between the 
components across the four fields is only ~70 kms~^ (see 
e.g. Fig. [2T]) . These values are on the order of the narrow- 
est line- widths measured, and are all much less than sound 
speed in the hot medium (^500 km s""*^ at ^ 10^ K), the 
escape velocity (-75-100 km s~^ iMartinll 19981 ). and the in- 
ferred velocity of the wind. This implies that turbulent ve- 
locities dominate over bulk motions in this whole 200 pc 
central region. 

We began this study with the aim of looking for ev- 
idence of the roots of the galactic wind, and to try and 
identif y where and how the flow is accelerated and colli- 
mated. [Heckman et all (|l995h and Martin (1998) identified 
and catalogued large-scale expanding supershells surround- 
ing the central regions of NGC 1569 that together form the 
large-scale galactic outflow. These authors find the velocity 
ellipses associated with shells A, B (on the western side) 
and E (on the north-eastern side) to have inner boundaries 
located at -10'', -30'' and -60" away from SSC A, re- 
spectively, placing the start of the collimated flows at least 
— 100 pc away from our IFU positions. We know that within 
the central 200 pc (as sampled by our observations), we do 
not see evidence for bulk motions consistent with the out- 
flow. This must therefore mean that the collimation process 
occurs between radii of 100 to 200 pc from SSC A. This 
interesting region clearly deserves further detailed study, as 
it must be very important in the galactic outflow formation 
process. 

Are we therefore detecting the galactic wind roots in 
our IFU positions? In IPaper ll we discussed how our obser- 
vations are reminiscent of the simulated gas motions within 
th e outer bounding sho c k of a n expanding zone as modelled 
by iTenorio-Tagle et"al] (|2006l ) , in the presence of a highly 
clumped medium. In light of this discussion and the re- 
sults presented above the answer must be yes, but in these 
inner regions, we are still well within the outflow's sonic 
point (equivalent in this case to the outer shock in the 
iTenorio-Tagle et aU clumpy H II region models) o r energy 
injection zone ([Shopbell k, Bland-HawthornI ll998l V and at 
this point the gas has yet to develop into an organised flow 
capable of reaching into the galactic halo. We therefore in- 
terpret the broad line as arising from gas that is mixing into 
the high sound- speed, hot gas outflow thereby mass- loading 
the wind. 

Finally we turn to Fig. [23l which shows the 



[Oiii]A5007/H/3 vs. [S ii]A6717+A6731/Ha nebular diagnos- 
tic diagram for all CI measurements in all IFU positions. On 
this plot, we als o show the theoretica l 'maximum starburst 
Une' derived bv iKewlev et al.l (|200lh . indicating a conser- 
vative theoretical limit for pure photoionization from the 
hardest starburst ionizing spectrum that can be produced. 
For points to exist above or to the right of this thresh- 
old, an additional contribution to the excitation from non- 
photoionization processes (e.g. shocks) is required. For a spe- 
cific metallicity or ionization parameter, however, the hard- 
ness of the ionization field is constrained and this limit may 
change, but since the model predictions are quite uncertain 
it is difficult to know by how much. For this reason we 
have employed the more robust 'maximu m starbu rst line' 
throughout this work, as is discussed in IPaper t Fig. [23] 
shows that spaxels with line ratios consistent with shocks 
are found only in positions 2 and 4. However, what is not 
immediately obvious from the maps show ing the spatial lo- 
cation of these points (fig 11 in IPaper B for position 2, and 
Fig. [18] for position 4) is that the distribution of their ra- 
tios are offset, with higher [S il]/Ha ratios found in position 
4. Since [S 11] emission can be strongly enhanced by shocks 
(|Dopital I1997I : lOev et all I2OOOI ). this position contains the 
strongest evidence for shock-excited emission. 

4.1 The evacuated cavity 

In our discussions of the individual fields, we found that the 
low surface-brightess inner cavity surrounding SSCs A and 
B can be identified at the edges of the FoV of positions 1, 
3 and 4, so we should be able to make a rough assessment 
of the conditions in this cavity by combining the data from 
these small sections. Firstly, we find that all the associated 
corners or sections emit only faintly in Ha, as would be 
expected from inspection of Fig. [1] Furthermore, all these 
areas have stationary or redshifted velocities in the range 0- 
30 kms~^ with respect to v^ys (with the exception of C3 in 
the south-western corner of position 1, which has redshifted 
velocities of up to 100 kms~^). This result is consistent with 
the evacuated nature of this region, implying that here we 
are seeing further into the galaxy to a region where the pre- 
dominant gas motions are directed away from us. 

In positions 1 and 4 on the north and west of SSC A, 
we observe very broad line widths in both CI and C2 for ar- 
eas associated with the cavity (up to 100 kms~^ in CI and 
250 kms~^ in C2), however the opposite is true in position 
3, where we find the narrowest CI line- widths of any field. 
Unfortunately, since our sample is very small and limited in 
spatial extent we cannot infer anything further. Additional 
observations are needed to better understand the gas condi- 
tions in this faint cavity. 



5 SUMMARY 

The small size of the CMOS IFU FoV (50 x 35 pc), cou- 
pled with the high spatial- and spectral-resolution of the 
instrument, have allowed us to examine each of the four 
fields observed in unprecedented detail, and relate the phys- 
ical conditions of each field to the general state of the ISM 
in this central region. We now summarise the main points 
arising from our study. 
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Figure 23. Dereddened CI flux ratios of [O ill] A5007/H^ vs. [S ll] A6717+A6731/Hq; for all four IFU positions. The dashed curve represents 
the theoretical maximum starburst line (|Kewlev et al.lEoOll ). 



• We see knots and filaments of gas which appear to be 
located at different lines-of-sight and expanding in different 
directions, attesting to the highly disturbed and fragmented 
nature of the ISM in NGC 1569. 

• We observe a narrow (^^35-100 kms~^) and broad 
('-^100-400 kms~^) component to the Ha line across all four 
regions. The total radial velocity spread between the com- 
ponents is only ^70 kms~^. 

• We also identify a third, high- velocity component in 
some regions, which we can confidently associate with shell 
expansion in positions 2 and 4, and po ssibly with the 
hi^h v elocity kinematic system identified bv lHeckman et aP 
([l99i). 

• We confirm an d expand on our conclusions of Westmo- 
quette et al. f 2007b. [PaperB ) for the origins of the two main 
components. The most likely explanation of the narrow com- 
ponent (CI) is that it represents the general disturbed opti- 
cally emitting ionized ISM, arising through a convolution of 
the stirring effects of the starburst and gravitational virial 
motions. This gives a region- wide characteristic turbulent 
width to CI (referred to as FWHMo), which is likely aug- 
mented by a varying contribution from unresolved kinemat- 
ical components (expanding shells, knots of gas) along the 
line-of-sight. FWHMo varies between 30-45 kms~^ across 
the four fields, and is in good agreement with the corre- 
sponding integrated H i velocity dispersions. 

• Our data support the conclusion of IPaper ll that the 
broad component (C2) results from the highly turbulent ve- 
locity field associated with the interaction of the hot phase 
of the ISM (the cluster winds) with cooler gas clumps - par- 
ticularly in position 1. In our model, the broadest emission 
is associated with the direct impact of the hot, fast-flowing 
cluster winds on t he cool gas knots, setting up turbulent 
mixin g layers (e.g. iBegelman k FabianI Il990l : [sI avin et al. 
Il993l ). Broad emission across the rest of the field most 



likely arises from material that is photo-evaporated or ther- 
mally evaporated from these layers through the action of 
the strong ambient radiation field pervading the central re- 
gions of NGC 1569, or mechanically ablated (stripped) by 
the impact of the winds. 

• We find surprisingly little evidence for shocked line ra- 
tios. This may be because the signatures of shocks are over- 
whelmed by the strong photoionization in the central re- 
gions. Positions 2 and 4 show the most convincing evidence 
for shocks, where in position 4 the emission is associated 
with a large expanding shell complex extending to the east 
of the FoV. 

HST imaging (Fig. [T]) clearly shows that the central re- 
gions of NGC 1569 are dominated by the effects of the cluster 
winds from SSCs A and B. This is exemplified by the swept- 
back, cometary-tail appearance of some of the remaining gas 
knots, and the existence of a large cavity surrounding the 
clusters, ^200 pc in diameter. 

This investigation has focussed on the small-scale char- 
acteristics of the ionized ISM projected within this central 
starburst zone, within which we find clear signatures of the 
multiphase nature of the ISM. We interpret the bright nar- 
row emission lines as arising on photoionized surfaces of 
clouds, and find them to be kinematically associated with 
the Hi. While some indications of normal expanding shells 
also are detected, these are relatively isolated and too slow 
to be important in driving the starburst wind. The ener- 
getic nature of this region reveals itself through the ubiqui- 
tous presence of broad components of optical emission lines 
with velocity FWHM ranging from 100 to ?^400 kms~\ The 
characteristic velocities of this material correspond to tem- 
peratures of 10^ - >10^ K, far higher than the -10^ K 
associated with the emission lines that we observed. 

The velocity offsets between broad and narrow compo- 
nents are modest, <50 kms~^, meaning that there is little 
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evidence for significant ordered gas flows within the central 
zone of the NGC 1569 starburst, where energy is being in- 
jected into the outflow. The flow dominated wind in NGC 
1569 must therefore form well beyond the region containing 
the massive star clusters. We interpret the combination of 
single peaked broad optical emission line components and 
the small net radial velocity offsets as arising from gas that 
is mixing into a high sound-speed, hot gas phase, whose 
presence dominates X-ray images of NGC 1569. The lack 
of bulk motions implies that this region is in a quasi-static 
state , as is expected for the inner regions of ga lactic winds 
fe.e:. iLeer Holzerlll980l : ISuchkov et al.l 1 19961 ) which pass 
through a sonic point. In this model the sonic point lies be- 
yond the regions we observed, at distances of > 100-200 pc 
from the starburst. In a forthcoming contribution (West- 
moquette et al. in prep.; Paper III) we will examine these 
conclusions in the context of the large-scale galactic outflow 
using deep Ha imaging and IFU observations. 

Finally, we consider the wider implications of our find- 
ings. Underlying broad components to Ha emission lines 
have been reported in a number of studies of giant extra- 
galactic Hii region s and violent s tarburst regions in galax- 



(e.g. IChu fcK ennicutt 19 ^: iMendez EstebanI Il997l : 



iHomeier k Gallagher.1999. : .Sidohet al.ll2006l ). Recently, we 
have added M82 to this list; from FST/STIS spectroscopy 
we find a ubiquitous broad Ha component throughout the 
starburst core Westmoquette et al. (2007b). The properties 
of these broad components resemble those we have found in 
NGC 1569 in that, compared to the narrow component, they 
have similar ionization conditions, radial velocities, and are 
spatially extended over the star -forming knots which con - 
tain young massive clusters (e.g. iMendez Estebanlll997l V 
It therefore seems likely that in these systems too, the broad 
emission line component is associated with the impact of 
cluster winds on cool gas knots. This suggests that the 
broad component may well serve as a powerful diagnostic of 
the galactic wind phenomenon since it probably traces both 
mass loading of wind material, as shown here, and mass en- 
trainment. The relationships between the properties of the 
broad component, the starburst regions, and the galactic 
wind flow merits further observational study. 
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